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Abstract 

This paper concerns stochastic perturbations of piecewise-smooth ODE systems relevant 
for vibro-impacting dynamics, where impact events constitute the primary source of ran¬ 
domness. Such systems are characterised by the existence of switching manifolds that divide 
the phase space into regions where the system is smooth. The initiation of impacts is cap¬ 
tured by a grazing bifurcation, at which a periodic orbit describing motion without impacts 
develops a tangential intersection with a switching manifold. Oscillatory dynamics near reg¬ 
ular grazing bifurcations are described by piecewise-smooth maps involving a square-root 
singularity, known as Nordmark maps. We consider three scenarios where coloured noise 
only affects impacting dynamics, and derive three two-dimensional stochastic Nordmark 
maps with the noise appearing in different nonlinear or multiplicative ways, depending on 
the source of the noise. Consequently the stochastic dynamics differs between the three 
noise sources, and is fundamentally different to that of a Nordmark map with additive 
noise. This critical dependence on the nature of the noise is illustrated with a prototypical 
one-degree-of-freedom impact oscillator. 


1 Introduction 

Many vibrating mechanical systems experience undesirable impacts that cause wear or sub- 
optimal performance. Occasional impacts may be permissible if they result from running com¬ 
ponents of the system at high speeds for greater efficiency, and some impacts are unavoidable 
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such as those due to random or rare events. In these cases it is important to have a clear un¬ 
derstanding of the dynamical behaviour that impacts may induce. Impacting dynamics is often 
complicated or chaotic because impacts are highly nonlinear phenomena [HE] GUIS IS]. 

For instance tubes in heat exchangers vibrate at high fluid velocities and impact against 
baffles in place to guide the fluid flow. Simple mathematical models of heat exchangers reveal 
that chaotic dynamics may be created at the onset of recurring impacts M- Rotating cutters 
spun at high speeds experience repeated contact loss with the material being cut. The resulting 
impacts between the cutter and the material may similarly induce chaotic dynamics PEI Sol. 
In contrast, some mechanical systems use impacts to achieve their function. Atomic force micro¬ 
scopes measure surface topography and the chemical properties of a sample on the nanoscale by 
gently hitting the sample with a vibrating cantilever. In this context it is important to under¬ 
stand the complex impacting dynamics so that the inverse problem of describing the sample can 
be performed effectively miss sg. 

Impacts can often be modelled accurately by carefully describing the deformations that com¬ 
ponents of the system undergo during impacts [HI US]- However, for the purposes of under¬ 
standing vibro-impacting dynamics, such a modelling approach is too cumbersome and a low 
degree-of-freedom ODE model can be more useful. Despite the low-dimensionality, such models 
have been shown to quantitatively match the experimental data of a variety of impacting systems. 
Examples include a cam-follower system involving occasional detachments between the cam and 
the follower mm, a pendulum experiencing near-instantaneous impacts with a solid wall [IS] , 
and compliant impacts of a steel block with an elastic beam nmm eh. 

Using a low degree-of-freedom ODE model, the evolution of the system between impacts is 
tracked in phase space. Periodic behaviour without impacts corresponds to a periodic orbit in 
phase space that does not reach the switching manifold, £, Fig. OJA. Physically, £ corresponds 
to locations or instances where mechanical components come into contact, or lose contact. As 
parameters vary, the system may transition from an impact-free regime to the repeated (though 
not necessarily regular) occurrence of impacts. In phase space, the transition occurs when the 
periodic orbit of the ODE model attains an intersection with £, Fig. [IJB. This is known as a 
grazing bifurcation. 

In this paper we study grazing bifurcations of the three-dimensional piecewise-smooth ODE 



Figure 1: Sketches of phase space illustrating a regular grazing bifurcation occurring at r/ = 0, 
where rj G M is a system parameter. 
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system, 



f L (u, v,w, rj) , 
f R (u,v,w,ri ) , 


u < 0 
u> 0’ 


( 1 . 1 ) 


where /l and f R are smooth functions, r) G R is a parameter, and the coordinates (u, v, w) are 
chosen so that £ is simply the coordinate plane u = 0. We assume that for rj < 0, there exists an 
attracting periodic orbit describing non-impacting dynamics located entirely in the region u < 0, 
and that the periodic orbit grazes £ at the origin when rj = 0. In the context of vibro-impacting 
systems, u < 0 corresponds to not-in-contact dynamics governed by /l, and u > 0 corresponds to 
in-contact dynamics governed by f R . Impacts may instead be modelled as instantaneous events 
with energy loss and velocity reversal, in which case a map is usually defined on the switching 
manifold to describe the action of an impact |22j . 

Theoretical studies of piecewise-smooth and hybrid dynamical systems have led to a useful 
classification of grazing bifurcations [23]. This paper concerns regular grazing bifurcations. The 
grazing bifurcation of (II. ip at rj = 0 is said to be regular if 


sgn(e]f L (0,v,w;ri)) = sgn (eJf R (0,v,w;i])) , 


( 1 . 2 ) 


for all (■ v,w,r /) in a neighbourhood of (0, 0; 0). This condition arises naturally in mechanical 
systems with compliant impacts, and implies that £ is neither attracting nor repelling at any 
point. 

As indicated in Fig. Q3-C, the steady-state dynamics of (11,11) for rj > 0 is often complicated. 
For this reason it is valuable to study the oscillatory dynamics using a return map based on the 
points on a Poincare section. A normal form for such a map for regular grazing in R 3 is the 
Nordmark map 


■i'i ■ l 

y %+1 


f g(x%,y %), 


| g(x uyi )- 

0 


Xi < 0 
Xi > 0 


(1.3) 
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with r, 5 G R, and x = ±1; as determined by the sign of certain coefficients, see [|2j The 
coordinates (x,y) represent points on a Poincare section of (11.11) . with fi G R the bifurcation 
parameter and the grazing bifurcation occurring at /./ = 0. The Nordmark map (II.3p . applicable 
also to models with instantaneous impacts [2T, 25] , includes only the leading order terms of the 
return map, so is valid only for dynamics close to the grazing bifurcation [23]. 

Each iteration of (jl.3p corresponds to one oscillation of (11.11) near the grazing periodic orbit. 
The utility of (1 1.3 p lies in the fact that the nature of the dynamics can be identified by the 
location of the corresponding points in the (/i, x)-plane. Specifically, (II.3p is formulated so that 
if Xi < 0 then the oscillation lies entirely in u < 0, and if Xi > 0 then the oscillation enters u > 0. 
The non-impacting, attracting periodic orbit shown in Fig. [TJ-A corresponds to the fixed point of 
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(1.5) 
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with x L < 0 and fi < 0. 

More generally, a periodic orbit of (II.ip appears as a finite set of points in the (fi, x ) plane, 
as in Fig. [2] which shows a typical bifurcation diagram of (j 1.3 jl . This figure shows a period¬ 
incrementing cascade, corresponding to several different periodic orbits, and apparently chaotic 
dynamics, as indicated by a cloud of points. For alternate values of r and 5, the fixed point 
(. x L ,y L ) may bifurcate directly to chaos [22] 26]. The square-root term in (II.3p is an artifact 
of the tangency between the periodic orbit and the switching manifold of (II.ip at the grazing 
bifurcation, and is responsible for the distinctive shape of the bifurcation diagram near fi = 0. For 
some vibro-impacting systems it is more appropriate for the return map to be piecewise-linear 
and either continuous 123 or discontinuous [28 j . Such maps predict fundamentally different 
bifurcation structures to those of (1 1 .31) . 

In order to properly explain complicated vibro-impacting dynamics, the effects of randomness 
and uncertainties needs to be taken into account. Mechanical systems are subject to background 
vibrations and other sources of noise. Experimentally measured parameters involve error, and 
some physical features are left unmodelled. For instance, one-degree-of-freedom models do not 
capture high frequency modes that are often excited by impacts [29]. 

To quantitatively describe stochastic impacting dynamics, stochastic averaging methods have 
proved useful for vibro-impacting systems that experience a wide range of impact velocities 
[3D, EU 021 Ell • If only low-velocity impacts are relevant, then it is useful to study (1 1.3 p . In 
his PhD thesis [03], Griffin studied (11.3p in the presence of additive white noise. He found that 
noise blurs bifurcation diagrams and washes out high-period solutions in the same manner as for 



h 

Figure 2: A bifurcation diagram of the Nordmark map (11.31) with r ps 0.5813, 5 ~ 0.1518, and 
y = 1. These parameter values correspond to the vibro-impacting system described in [J3], with 
(fcosc, b osc , fcsupp, frsupp, d) = (4.5, 0.3,10, 0, 0.1). The fixed point of (JT3]) for /i < 0, given by dH5]), 
corresponds to an attracting, non-impacting periodic orbit of period 27T. The map (11.3p has 
attracting 3 and 4-cycles for certain values of /r > 0 as shown, as well as an attracting 5-cycle in 
the approximate range, —0.00004 < fi < 0.00041. These correspond to periodic orbits of period 
approximately equal to 2/c7r, for k = 3,4,5, that experience one impact per period. There also 
appears to be a chaotic attractor for the approximate range, 0.0135 < p < 0.0155. 


4 

















smooth maps, such as the logistic map [35, 36]. Recently it was shown that white noise added to 
(II.ip translates to additive white noise in (11.31) . [ 37] , Such a noise formulation may be sensible for 
vibro-impacting systems for which a forcing term or external fluctuations represent a significant 
source of uncertainty. 

However, impact events themselves constitute a substantial source of randomness. The pur¬ 
pose of this paper is to construct and analyse stochastic versions of (1 1.3 [) for which randomness 
stems purely from impact events. We consider three different types of impact noise for (11.11) 


using an Ornstein-Uhlcnbeck process with stationary density J\f 


0 — 

2 v 


, where e <C 1 represents 


the noise amplitude and v > 0 is the correlation time. We first consider uncertainties in £, then 
uncertainties in fn, and lastly uncertainties in fj .j in the white noise (zero correlation time) limit. 

Here let us indicate the forms of the stochastic maps that we obtain. Coloured noise in £ 
leads to random perturbations in both the map and the switching condition as 


%i -\-1 

Vi +1 


Ni(xi,yi) 


{g{xi,yi) , 


I g(xi,yi) - 

Xy/xi + «i 6 

0 


Xi + Ki & < 0 


Xi + Ki£i > 0 


( 1 . 6 ) 


where are Gaussian random variables, and > 0 is a constant. For coloured noise in the 

impacting dynamics we obtain 




'g{xi,yi ), 


Xi -\-1 
Vi+ 1. 

= N 2 {xi,yi ) = < 

g(xi,yi) - 

X ^2 (ti)y/xl 

0 


(1.7) 


for a particular nonlinear function k 2 . For white noise (v —> 0) in the impacting dynamics the 
map takes the form 


%i -\-1 

Vi +1 


N 3 (xi,yi) 


{9(xi,Vi) 

g{xi,yi ) 


x«3 (ri,hi)y/xi 

0 


+ Ki[hi)Xi , 


Xi < 0 
Xi > 0 


( 1 . 8 ) 


where r; and hi are random variables, and k 3 : M 2 —> R and : M —> M 2 are nonlinear functions. 
Notice that N± is stochastic for Xi < 0, whereas N- 2 and N 3 are not. This is because noise in 
£ generates anomalous crossings of Xi = 0. Consequently N\ exhibits stochastic dynamics for 
fji < 0, while N 2 and N 3 do not. N 2 and N 3 involve noise terms proportional to y/xi, and for this 
reason exhibit increasing variability for larger values of /i > 0. 

To obtain a more detailed comparison of Ni, N 2 and N 3 , we first carefully derive (11.3ft in )|2j In 
(J3]we introduce a prototypical compliant vibro-impacting system to illustrate our results. In )|4] 
we add randomness and derive (1 1. 61) - (1 1. 81) . Since (11.61) - (1 1. 81) involve fundamentally different noise 
terms, they exhibit different sensitivities to the noise amplitude e. Therefore we use different 
values of e in each of the different models, in order to make appropriate comparisons. For each 
Nj (j = 1,2,3) we write £ = e^a and identify the appropriate value £j, so that N 3 , N 2 and N 3 
display roughly the same dynamics for the vibro-impacting system with the parameters of Fig. [2] 
and y, = 0.03 (chosen for illustration) and a = 1. This enables us to quantitatively compare , 
N 2 and N 3 in Conclusions are presented in §|6l 
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2 A derivation of the Nordmark map 

In this section we derive a return map for the generic deterministic system (II.ID valid near the 
grazing bifurcation rj = 0 and provide an explicit coordinate change that transforms the map to 
(II.3p . Such a derivation is given in [23] [38]. We provide an explicit construction of (II.3(1 here in 
order to provide a basis for deriving the stochastic maps in §j4j Throughout this section we write 
z = ( u , v , w) for points of (II.ip in R 3 . 

As discussed in ]JT] we assume that with r/ = 0 the system (11.11) has a periodic orbit that 
intersects the origin, but is otherwise contained in u < 0. This implies that u = 0 for fi at the 
origin with 77 = 0 (i.e. e|/b(0, 0, 0; 0) = 0). For simplicity, we assume that we further have 

eJf L (u,0,w\ri) = 0, (2.1) 

for all (u, w ; r/) in a neighbourhood of (0, 0; 0), which can usually be imposed by an appropriate 
coordinate change. This assumption is particularly useful in the case that randomness is present 
in the switching condition, see 1 14.11 

The key to deriving the Nordmark map (II.3D is selecting a convenient Poincare section and 
constructing a discontinuity map that accounts for the difference between fi and fn, that is, 
the difference between impacting and non-impacting dynamics. We let II represent a generic 
Poincare section of (II.ip that lies in u < 0 and intersects the grazing periodic orbit transversally, 
and let II' be the coordinate plane v = 0, see Fig. [3] 

Given zq G II, we define z 1 G II' as the next intersection of the forward orbit governed by 
fi (i.e. ignoring the switching condition at u — 0) with II'. If u > 0 for the point z±, then z\ 
does not represent the true intersection of the orbit of (II.ip with II. Nevertheless, we study the 



Figure 3: A schematic diagram illustrating the construction of the Nordmark map (11.31) for the 
three-dimensional piecewise-smooth system (11.11) . The solid curve represents an orbit of (II.ip . 
The dashed curves show virtual extensions of this orbit into u > 0 as governed by jj J . 
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return map on II' (z\ —> zp) rather than the return map on II (zq —> Z 5 ), because the two maps 
are conjugate and the map on II' has a simpler form than the map on II. 

To derive the return map on IT, we must consider three additional points, z 2 , z 3 , and 24 , see 
Fig. [31 The points z 2 and z 3 correspond to the entry and exit locations of the orbit with the 
impacting region, u > 0. Then z\ is obtained by travelling forward from the entry point z 2 to II' 
using f L , whilst 24 is obtained by travelling backward from the exit point £3 to II' using f L . The 
discontinuity map is defined as 24 = D(z\). If u < 0 for the point Zi, then 24 = Z\. 

In order to derive an explicit expression for the discontinuity map, we consider the three 
steps, Z\ —y z 2 , z 2 — y z 3 , and z 3 —> Z4, individually, and expand /l and fn about r/ = 0 and the 
origin by writing 


fj(u,v,w-,rj) 


ajv + 0(2) 

~Pj + 0(1) 
7 J + 0(1) 


( 2 . 2 ) 


where aj,Pj,^/j G M, and J = L,R. I 11 02.21) . and throughout this section, O(k) is short-hand 
big-0 notation for O (ji/u, v, w, r]) k ^j . Notice that yfu is assumed to be of the same order as v, 
w, and 77 , which is appropriate in view of the tangency between the grazing periodic orbit and 
£. We assume 


a L, OIr, Pl, Pr > 0 , 


(2.3) 


such that the tangency is quadratic for both f L and _/Vi, and has the orientation depicted in 

Fig. El 

From series expansions of the orbits governed by /l and fa, see [23108], we obtain the following 
formulas for the three steps in the discontinuity map in the case U\ > 0 (writing z t = (ut, v il Wi)). 
For Z\ y z 2 . 


for z 2 -> z 3 , 


and for z 3 —> 24 , 


V2 


W 2 



VWE 




y/U\ + 0(2) 


Wi - 


VZ'YL 

V a L^L 


sju{ + 0(2) 


V3 


_W 3 _ 


U4 


W4 



—v 2 + 0(2) 

W 2 + ^v 2 + 0(2) 


+ 0(3) 


W 3 + J^v 3 + 0(2) 

By combining (I2.4MI2.6I) we get (« 4 , u^) in terms of (ui,Wi), 


(2.4) 


(2.5) 


( 2 , 6 ) 


where 


U4 


U\ + 0(3) 

W4 


wi - Cy/TR + 0(2) 


2y/2P L fjL__ 7 iA 

\/&L \Pl PrJ 


(2.7) 


( 2 . 8 ) 
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The discontinuity map is then 


D(u, w; rj) 


u 


w 


u + 0( 3) 

W — Cyfu + 0(2) 


u < 0 


u > 0 


(2.9) 


To complete the map (11.31) . which represents Z\ —> z e as shown in Fig. [3j we must combine 
D(z\) with the global return map z 6 = G(z 4 ). G depends on global properties of /l and is 
smooth, so we can write 


for some 


G(u, w, rj) 



+ br/ + 0 (( u,w,r]) 2 ) , 


A = 

an 

«12 

, b = 

'bj 


«21 

022 


b 2 


( 2 . 10 ) 


( 2 . 11 ) 


where each a % j , hi G M. Then the desired return map on II' is the composition G o D. 

Finally we apply a coordinate change to convert the map to the normal form (j 1.3 [) that 
involves only three parameters, r, 5 and x- Under 


X 

l 

1 0 

0 


U 

y 

„ 2 „2 

— 022 0-12 

bi 


W 

y 

u 12 c 

0 0 

(1 — 022)^1 + ®12^2 


y 


( 2 . 12 ) 


and with higher order terms omitted, G o D transforms to (II.3p with 


r = trace(/l) , 6 = det(A) , x = s g n ( a i 2 c) . (2.13) 


Since the Nordmark map keeps only leading order terms for (u, w\ r]) near (0, 0; 0), terms that are 
linear in x are omitted since they are dominated by y/x : refer to |39] for a further discussion. In 
(I2.12p we require a 42 ^ 0 and c ^ 0. These represent non-degeneracy conditions for the grazing 
bifurcation. 


3 An oscillator with compliant impacts 

To motivate and illustrate our results for stochastic Nordmark maps, we consider the prototypical 
vibro-impacting system shown in Fig. [4] and studied in [371 H0l ITT| . This system consists of a 
harmonically forced one-degree-of-freedom linear oscillator that experiences compliant (or soft) 
impacts with a support, and we use the following non-dimensionalised equations to model the 
dynamics: 


—k OS c( u + 1) — b osc ii + Fcos(t) , u < 0 

^osc(^ + 1 ) (bosc + frsu PP )h ^supp(^ T dj T F cos (t) , u 0 


u = 


(3.1) 



























Here u(t ) denotes the location of the block, which has the equilibrium position u = —1. A 
rigid stop prevents the support reaching a position with u < 0, and prestresses the support by a 
distance, d > 0. The constants k osc , b osc , k supp and b supp represent non-dimensionalised spring and 
damping coefficients for the oscillator and support. We neglect the mass of the support, ignore 
energy loss at impacts, and assume that whenever the block is not in contact with the support, 
the support is located at u — 0. Experiments of simple vibro-impacting systems with compliant 
impacts have shown that low-dimensional models such as (13. ip can quantitatively match the 
physically observed dynamics near grazing bifurcations mmm- 

Here we treat the forcing amplitude, F > 0, as the primary bifurcation parameter, and assume 
kosc, b osc > 0. The steady-state solution (behaviour in the limit t —> oo) to (13.IK with u < 0 is 


U ss{t) 


, {ko sc - 1) cos (t) + b osc sin (t) 

(k osc - l) 2 + 6L 


(3.2) 


When F < _F graz , where 

Eg r a Z = \J(k osc — l) 2 + 6 2 sc , (3.3) 

the maximum value of u ss (t) over one period is negative, and so u ss (t) is an attracting non¬ 
impacting periodic orbit of (13.ip . The critical value, F = F graz , is a grazing bifurcation at which 
u ss (t) has unit amplitude and attains the value u = 0 at times t = £ graz + 27r/c, for fceZ, where 



and 0 < t graz < 7r. 

To convert (13.IK to the general form (II. IK . we dehne 

v = it , w = (t mod 27 t) - t graz , r) = F - F graz . 


(3.4) 


(3.5) 


Here the phase space of (13.IK with (13.51) is isomorphic to M 2 x T, rather than M 3 , but this does 
not affect the bifurcation structure near grazing. For (13.IK with (13. 5K . the coefficients in (12.2K . 


stop 



d 


u = —1 u = —d u = 0 

Figure 4: A schematic diagram of the vibro-impacting system modelled by (13. ip . This system 
exhibits a regular grazing bifurcation with a square-root singularity because with k supp > 0 and 
d > 0 the equations of motion are discontinuous at the grazing point. 
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which describe the behaviour of the system near the grazing point, are given by 


®L = 1 , @L = 1 , lL = 1 , 

1 5 /3r 1 “h ^supp^ j 1 5 


(3,6) 


and by evaluating (|2.8|) with (|3.6|) we obtain 


c = 


2 \/2 A^supp d 


1 + &supp^ 


(3,7) 


In addition, from the general solution to (13. 1 [) we find that the coefficients in the global map 
(12. 10p are given by 


(3,8) 


f 2 tt 

0 

1 

) , b- 1 

1 - a n 

V 

^OSC 

^osc 

) F 

/ graz 

— a 21 


4 Incorporating randomness into the Nordmark map 

To model noise and uncertainties we use the one-dimensional Ornstein-Uhlenbeck process 

d£(t) = -~£(t) dt + — dW(t ) , (4.1) 

v v 

where e, v > 0 are constants and W(t) is a standard Brownian motion. Given an initial value 
£(0) = f 0 , at any positive time £(£) is a Gaussian random variable with mean and variance 

E [^(^)l^(0) = ^o] = Co , Var[f(t)|f(0) = f 0 ] = — (l -e"^) . (4.2) 


In the limit t —> oo, £(t) ~ J\f 


0, 


£ _ 

2u 


, where Acr 2 ] denotes the Gaussian distribution of mean 


fi and variance a 2 . The correlation time of (14.11) . defined as / 0 °° with £(0) ~ M 


0 — 

2v 


is equal to v. 

In our context, £(£) is coloured noise and the parameter e governs the size of the noise. 
Unlike white noise, £(£) has an inherent time-scale, u, and is suitable for modelling various types 
of uncertainties in mechanical systems, such as background vibrations J32]. In the white noise 
limit, forcing by £(t) becomes a diffusion process edW(t). 


4.1 Stochastic switching 


We first consider the following stochastic perturbation of (II.ip . 



f L (u, v, w, rj) , 
f R (u,v,w,rj) , 


u + £,(t) < 0 

u + £(t) > o 


(4.3) 
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where £(i) is given by (14.ip . In (I4.3[i randomness is present in the switching condition, while 
evolution between switching events remains deterministic. We expect (14.3p to be applicable to 
a wide variety of piecewise-smooth systems. For the vibro-impacting system of £(f) may 
capture uncertainties in the point at which contact between the block and support occurs or is 
lost. For switched control systems, £(t) may correspond to measurement errors that produce 
variability in evaluations of switching rules (43, [44]. 

Here we consider orbits of (14. 3 p that are close to the grazing periodic orbit of (11.111 . Orbits 
of (14. 3 p near grazing only spend short lengths of time in the region u > 0 while passing near the 
origin, and for simplicity we suppose that the value of v is large compared to such times. In this 
case it is reasonable to approximate £(t) by a constant while an orbit is near the origin. With 
this approximation, the sum u{t) + £(£) does not switch sign more than twice as the orbit passes 
near the origin, which substantially simplifies our calculations below. 

We let denote the value of £( t ) during the z th instance that the orbit of (14. 3 p passes near the 
origin. The time between between consecutive traversals near the origin is well-approximated by 
the period of the grazing periodic orbit, call it T. With this approximation, 







(4.4) 


To derive the stochastic version of (11.311 for (14.3p with (14.411 . we first derive the induced 
stochastic discontinuity map. Here condition (12.ip is useful, as it implies that an orbit governed 
by fi attains a local maximum value of u at an intersection with n'. For u\ + £j < 0, we conclude 
that u{t) + & < 0 as the orbit passes near the origin, so that it 4 = u\. If instead U\ + G > 0, 
then the discontinuity map D is given by (12.7p except that u + £, : appears inside the square root 
because this quantity represents the distance from the switching condition. That is, 

W 4 

U>4 


w l 


U\ + 0(3) 

CyJ U\ + + 0(2) 


(4.5) 


where 0(k) = 0 ( (y/u, y/\£i\,v ,w,r )j I. By combining (14.5ft with the global map G, applying 


the coordinate change (12.12p . and dropping higher order terms, we obtain 


^ 2+1 
Vi +1 


Ni(xi, t/i) = 



r 

1 




0 



r 

1 


< 

-5 

0 



Xi 


+ 


0 

1 

Xi 


h , 


Xi + 


< 0 


y% - X\ Xi + 


+ 


h , Xi + 2 > 0 


(4.6) 


N\ is the stochastic Nordmark map corresponding to (14.3p . Notice that randomness in the 
switching condition of (14.3p has translated to randomness in both the switching condition of 
(11.31) and in the image of the map with Xi > 0. In contrast, a piecewise-linear map for which 
randomness is present purely in the switching condition is studied in [45]. 

In order to fairly compare N\ with other stochastic versions of (11.31) in [|5j we estimate the 
effective size of the stochastic contribution for our illustrative parameters values of the prototyp¬ 
ical system (13.11) and a representative value of /i = 0.03. This motivates us to express £ in terms 
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of a scaled parameter e = eqa, and to obtain comparable stochastic contributions for fixed a in 
the different cases. The square-root term of N\ is yjxi + where fiq = -4~n. If + is large 

a 12 c 

relative to K \£*, then this term is well-approximated by y/xl+ 2 ^ and so 2a 2 estimates 

the additive stochastic contribution to Ah- For the impact oscillator with the parameter values 
of Fig. [2] (here c = \/2 and au ~ 0.1227) and using Xi = 0.025 (corresponding to the value of 
Xi > 0 in Fig. [ 2 ] for fi = 0.03), this quantity is approximately 100£j. With v = 0.5 (used in 
§5]), the standard derivation of the stochastic contribution is approximately lOOe. Therefore, for 
e = £ia, where 

£i = 0.0001 , (4.7) 

the standard deviation of the stochastic contribution is approximately 0.01 when a — 1. 

To illustrate the accuracy of Ah, Fig. 0 compares iterates of Ah (black dots) with a numerical 
solution to (14.31) for the vibro-impacting system of 1|3] (purple dots) using £ = i 1 . For the given 
parameter values, the system has an attracting 3-cycle in the absence of noise. For this reason, 
both sets of points are grouped about the 3-cycle. The two sets of points are slightly separated. 
This is because the form of the deterministic Nordmark map does not include higher order terms 
of the true return map, as observed by the separation of the values taken by the deterministic 
3-cycles shown in Fig. [5l The size and shape of the spread of the two sets of randomly generated 
points are similar, as is their location relative to the deterministic values of the map. This 
demonstrates that Ah can accurately capture the stochastic dynamics of (14.31) . A more precise 
characterisation of the accuracy of Ah is beyond the scope of this paper. 
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Figure 5: A phase portrait comparing the stochastic return map Ni, (14.6ft . with a numeri¬ 
cal solution to (14.31) . The three groups of purple dots were obtained by numerically solving 
(14.31) with (14. If) . v — 0.5 and £ = ey, for the vibro-impacting system (13.lj) with (13.51) and 
(k osc ,b osc ,k supp ,b supp ,d) = (4.5,0.3,10,0,0.1) (as in Fig. EJ) and F F graz + 0.005558 (which 
corresponds to /i = 0.03). More precisely, 1000 points on IT (labelled Z\ in Fig. |3|) were obtained 
by numerically solving (14.3p . and these were transformed to (x, ^-coordinates by applying (I2.12p 
and (13.51) to produce the purple dots. The three groups of black dots are 1000 iterates of N { with 
(14.4p and parameter values chosen to match the vibro-impacting system (specifically, v = 0.5, 
£ — § 1 , T — 27r, 11 = 0.03, r ~ 0.5813, 5 ~ 0.1518, + = 1, c = \/2 and au ~ 0.1227). The 
deterministic 3-cycles of (14. 3 p and (14. 6 p are shown with triangles. 
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4.2 Additive coloured noise with a large correlation time 

Next we consider the case where randomness and uncertainty in (II.ip is associated with /#. For 
mechanical systems with impacts, this corresponds to variability in the evolution of the system 
during an impact. For simplicity we include noise in only the u-component of fn, that is 




' f L (u,v,w,rj) , 


u 

V 

w 

= < 

f r(u, V, W\ Tj) + 

1 - 

o S 

-Uy> 



< 

- 1 

O 


(4.8) 


where £(t) is given by (14.ip . Indeed, for the vibro-impacting system of Fig. [4J if noise is incor¬ 
porated into the force on the block when it is in contact with the support, then the equations of 
motion may be put in the form (14.81) . With noise added to the w-component of /r (or A), orbits 
may cross u = 0 many times in a short time frame which makes the system substantially more 
difficult to analyse. We leave such considerations for future work. 

As in 1 14.11 we consider near-grazing orbits and assume that the value of v is much larger 
than the time each orbit spends in the region u > 0. In this case it is reasonable to treat 
£(t) as constant while u > 0. During the z th instance that an orbit passes near the origin, we 
approximate £( t ) by £j, distributed according to (I4.4j) . In this scenario the three components of 
the discontinuity map (12.4ft - (12.6p are unchanged except that /3r is replaced by /3r — in (12.51) 
(because the u-component of the system with u > 0 is given by — (3r + & + 0(1), see (I2.2p h By 
combining (12.4h - fl2.6h we fold that for u\ > 0 the discontinuity map is given by 


W 4 

U>4 


U\ + 0(3) 


W\ — c 


TL 1R. \ 

)V^+o( 2) 

P L P R J 


(4.9) 


and therefore the corresponding stochastic Nordmark map is 


%i -\-1 

Vi+1 


= N 2 (xi,yi) = < 


r 1 
-h 0 


T 


+ 


h ; 


Xi 


Vi - X 


JR 


T L _ 

Pl Pr-U 

il 'm ~ 

Pl Pr 


+ 




Xj < 0 


Xi > 0 


(4.10) 


TL . 


JR 


With 

2 ^ 
2-Si ~ 


Notice that with ^ = 0, N 2 is identical to (II,3jl . 

We can write the stochastic component of N 2 as x« 2 (£i)y/xi, where K 2 (£i) = 

Pl Pr 

the parameter values of the impact oscillator (13.6ft . and k supp d = 1, we have k 2 (^) = 2 — ^ 

1 — y- Therefore the noise provides a multiplicative stochastic contribution of approximately 
ignoring signs. In order to compare the effect of the noise to the other cases, we write 
£ = e 2 a, choosing e 2 so that the standard deviation of the stochastic contribution is 0.01 when 
a — 1. For Xi = 0.025 and v = 0.5 (as in 1 14.1ft . the standard deviation of j s approximately 
0.08e, so we therefore choose 

e 2 = 0.125 . (4.11) 
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Fig. [ 6 ] compares iterates of N 2 to intersections with II' of a numerical solution to (14. 8 p using 
e = § 2 - As expected the two sets of points are similarly distributed about the deterministic 
3-cycle. 


4.3 Additive coloured noise with a small correlation time 


Lastly we consider (I4.8ji in the white noise limit, v = 0. In this case (I4.8ji reduces to a diffusion 
process forced by white noise, specifically £(t) dt is replaced by e dW(t). By using (12.2(1 to expand 
fa, (14.8p for u > 0 may be written as the three-dimensional stochastic differential equation 


du(t ) 
dv(t ) 

— 

a R v(t) + 0 ( 2 ) 
-fin + 0 ( 1 ) 

dt + £ 

"o' 

1 

dW(t) . 

(4.12) 

dw(t) 


7 r + 0 ( 1 ) 


0 




To obtain a stochastic Nordmark map corresponding to this scenario, we first derive the stochas¬ 
tic version of the middle component of the discontinuity map (12.511 . Given an initial point 
(u(0), u(0), ic(0)) = (0,v 2 ,w 2 ), where v 2 > 0 and w 2 E K are small, the desired values of n 3 and 
W 3 are given by the point (Cfi^uq) of first return for the stochastic process (j4.12(1 to u — 0. 
First return or first passage problems are an important class of theoretical problems in stochastic 
calculus with applications traditionally in finance and chemical kinetics [4Bl, W71 148] . 

We approximate (j4. 12 jl by keeping only the leading order contributions, i.e. 


du(t) 
dv(t) 


a R v(t) 

—fin 


dt + e 


dW(t) , 


«(0)' 


" 0 " 

."(0). 


? 2 _ 


(4.13) 


together with w(t) = w 2 + 7 Rt . With this approximation we are able to provide an explicit 
expression for the joint probability density function of the return location and time. A formal 
justification for the omission of the higher order terms is left for future work. 
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Figure 6 : A phase portrait comparing the stochastic return map N 2 with a numerical solution 
to (14.8p . The three groups of purple dots were obtained by numerically solving (14. 8 p with (14.ip . 
v = 0.5 and e = e 2 , for the vibro-impacting system (13.lh with (13.5ft using the same parameter 
values as in Fig. [5] The three groups of black dots are 1000 iterates of N 2 with 11.11) and 
parameter values matching those of the vibro-impacting system (refer to the caption of Fig. [5]). 
The deterministic 3-cycles of (14. 3 p and (14.6p are shown with triangles. 
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We introduce the change of variables 

Pr 


V = 


ol r v\ 


u 


Q = 


V 2 


Pr 

s = —t 
V2 


This puts (I4.13P in a standard form studied in 


dp(s) 

dq(s) 


q{s) 


ds + y/g 


dW(s) 


p(0)- 


"O' 

9(0). 


1 


where 


e = 


/3rv 2 


(4.14) 


(4.15) 


(4.16) 


Specifically, q is a diffusion process with constant drift, and p(s) = f* q(s) ds. Therefore p(s) 
may be interpreted as integrated Brownian motion with constant drift. 

We let J-(r, h] g) denote the joint probability density function for the first return of (I4.15|) 
to p = 0, at a time s = r > 0, and location q = —h < 0. In [5Uj, McKean derived an explicit 
expression for T in the case of zero drift by computing the inverse Kontorovich-Lebedev transform 
of the corresponding renewal equation. In |49] , Atkinson and Clifford extended this result to the 
case of non-zero drift by applying the Radon-Nikodyn derivative. Specifically 

(2^ I (r “ 2)2 “ 2 < r - 2 >( ft - !) + - l) 2 ]) ' ^ 


F(r. h\ i>) — exp 


erf 


ngr* 




(4.17) 


where erf(-) is the error function. The constant g > 0 governs the shape of T. The limit g —> 0 
corresponds to the deterministic case, for which r = 2 and h — 1. With a small value of g, T is 
roughly Gaussian. I 11 contrast, the limit g —* 00 corresponds to the case of no drift, as in [50], or 
to the limit v 2 —> 0. In this limit the marginal probability density function for r is asymptotically 
proportional to r~ 4 , for large r, and so is long-tailed [49] , 

In view of the scaling (14.1411 . the stochastic version of (12.5p corresponding to (14. 12)1 is given 
by 

V 3 
w 3 


—hiV 2 + 0(2) 

w 2 + ^v 2 + 0(2) 


(4.18) 


where and hi have the joint probability density function (14. 1 7[) . By combining (14.1811 with 
(12.4ft and (12.611 we obtain 


r n 


h'ju x + 0( 3) 

U4 



W4 


Wl-c( 2P i L 1™* ) y/u{ + 0 (2) 



V Pl p r J J 


(4.19) 


which represents the stochastic version of the discontinuity map for points with U\ > 0. Then 
using (14.19p we arrive at the following stochastic Nordmark map 


^2+1 
y *+1 


N 3 (xi,yi) = < 


T 1 

-5 0 


Xi 


Vi 


+ 


T + au(h] — l) 1 
-5hf 0 


h , 


Vi-X 


Xi 

( TL( fe i + 1 ) 7 B. r i \ 

WL R \ 

XL _ TR I 

P L PR J 


Xi < 0 

/i , Xi > 0 

(4.20) 
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In (I4.20p we treat each pair (r 8 , hi) as a two-dimensional stochastic random variable with proba¬ 
bility density function T(ri,hi\ Qi), where 


_ g 2 y^Z 
Ql /?r\/2/3l |a-i2c |\fx ’ 


(4.21) 


which results from combining (14. 16[) with v 2 ~ , (E3J), and ui = a 2 2 c 2 x, (12. 12)> . 

We now estimate the size of the stochastic contribution in N 3 . The leading order stochastic 

+ 1 ) 7 R r j 

component of N 3 is X K 3( r i,K)\/xi, where K 3 (r^hi) = — 2 '%_- La 2,a - R . With (13.6p and k supp d = 1, 

Pl Pr 

we can write k 3 (r*, hi)y/xl = y/xl + (hi — 1 — n j^ L ) y/xi, where the second term represents the 
multiplicative stochastic contribution, ignoring the sign of y, because the deterministic values of 
hi and r* are 1 and 2 respectively. 

With a small value of g (g < 0.03 is suitable), T(r, h] g) is approximately Gaussian because 
the effective noise amplitude in (14.151) is small. By (14.2ip . this approximation is valid when, 
roughly speaking, £ is not too large and Xi is not too small. From (14.17P we determine the 
covariance matrix of the Gaussian approximation to be 


Cov(r, h; g) 


2 q |"4 1" 

y i i ’ 


(4.22) 



and it follows that in this approximation the linear combination h— | has standard deviation 

The standard deviation of the stochastic contribution in N 3 is therefore approximately y^-y/Xi- 
Following the previous cases, we write £ = e 3 a, and choose e 3 so that the standard deviation of 
the stochastic contribution is 0.01 when a = 1. For N 3 , this quantity is — xfe ~ 0.46£, 

using (I4.2ip . the parameter values from Fig. [21 and Xi = 0.025. Therefore we let 


£ 3 = 0.022 . (4.23) 

Here gt ~ 0.008 when a = 1, and so for these values the Gaussian approximation to (14.171) is 
justified. 

Fig. [7] compares iterates of N 3 to intersections with n' of a numerical solution to (14. 8 p using 
£ = £3. As with the previous two figures, this shows that N 3 can accurately capture the stochastic 
dynamics of (I4.8p . 


5 Stochastic dynamics 

In this section we explore the dynamics of the three stochastic Nordmark maps, N \, N 2 and N 3 , 
and discuss how the different forms of these maps is evident in their dynamical behaviour. To 
briefly summarise, N 3 applies to the system with stochastic switching (14.31) . whereas N 2 and N 3 
apply to (14.81) . For N 3 and N 2 it is assumed that the value of u (the correlation time of the noise) 
is large relative to the times that orbits spend in u > 0. As these times are rarely larger than 
t — 0.05 for the parameter values considered here, we take v = 0.5 in N\ and N 2 to ensure that 
the correlation time is large enough. The values of G in N 3 and N 2 are distributed according to 
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Figure 7: A phase portrait comparing the stochastic return map N 3 with a numerical solution 
to (14. 8[) . The three groups of purple dots were obtained by numerically solving (14. 8 ft with (14.IF 
v = 0 (in which case £(t) dt is replaced by edW(t)) and e = £ 3 , for the vibro-impacting system 
(13.ip with (13. 5p using the same parameter values as in Fig. [5j The three groups of black dots 
are 1000 iterates of N 3 with (14. 1 7jl and (14.21)1 and parameter values matching those of the vibro- 
impacting system. The deterministic 3-cycles of (14.3p and (14.6(1 are shown with triangles. 


(14.4(1 . N 3 corresponds to the limit u —> 0. In N 3 , Ti and hi are distributed according to (14.17(1 
and depend on the value of Qi (14. 21 jl . 

For each Nj we have written 

e = £ja , (5.1) 

where the £j are given by ((4.7(1 . (14. lip and (14.23ft . These values have been chosen such that for 
a given value of a, the size of the stochastic contribution in Ni, N 2 and N 3 is roughly the same, 
at least when /i = 0.03. For n close to 0, the stochastic contributions are noticeably different for 
these choices of £ r 

In §5.1l we look at stochastic bifurcation diagrams in order to obtain a basic understanding of 
how the stochastic dynamics differs with the value of //. In the subsequent parts of this section 
we study two-dimensional invariant densities in order to gain a deeper understanding of the 
dynamics. 

5.1 The dependence of /i on the size of noise response 

Fig.[ 8 ]shows stochastic versions of the bifurcation diagram shown in Fig.[2]for the three stochastic 
Nordmark maps. As expected, the noise blurs the bifurcation diagram. In panel A, which 
corresponds to the map A r 1 , for values of /i very close to zero (say \n\ < 0.002) the points are 
relatively highly spread. This is because here the deterministic map has an attractor near x — 0, 
so in Ni the sign of Xi is often different to the sign of Xi + That is, the choice of the 

half-map of N 3 is regularly determined by rather than x % . Furthermore, as shown in 1 14.11 the 
leading order component of the stochastic contribution to the right half-map of N\ is inversely 
proportional to y/xi, which for very small /i is large relative to its value for fi away from zero. In 
contrast, with 0.03 < < 0.05 say, the underlying attracting 3-cycle is sufficiently far from x = 0 

so that the sign of x t rarely differs from that of Xi + The points are randomly distributed 
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Figure 8: Bifurcation diagrams of the stochastic Nordmark maps Ni (panel A), N 2 (panel B), 
and A 3 (panel C), with r ~ 0.5813, 5 ~ 0.1518, and x — 1- These parameter values correspond to 
the vibro-impacting system (13. l|h with (k osc ,b osc ,k suppi b supp ,d) = (4.5,0.3,10,0,0.1). The black 
dots are iterates of Ai, N 2 and N :i with transient points omitted. The noise amplitudes are given 
by (14.7)1 . (14. lip and (I4.23jl with a — 1, and v = 0.5 for panels A and B (panel C corresponds to 
v = 0). In each panel the deterministic bifurcation diagram (Fig. [2]) is shown in blue. 
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about the 3-cycle due to noise in the right half-map of N\. The size of the deviation decreases 
with increasing /r, because the strength of attraction of the 3-cycle increases with yu. 

In panels B and C, which correspond to N 2 and N 3 respectively, the bifurcation diagrams 
show no variability for yu < 0. This is because for /i < 0, ( x L , y L ) (jl .5jl is a fixed point of N 2 and 
N 3 . For 0.002 < /i < 0.01, the size of deviations about the underlying attracting 4-cycle increases 
with /i. This is primarily because the coefficient of the i/x-term of N 2 and N 3 is random, and 
for the 4-cycle this value of x increases with /i. For 0.03 < /i < 0.05, the size of deviations varies 
little with yu because the increased variability caused by a larger value of x is balanced by the fact 
that the strength of attraction of the 3-cycle increases. Panels B and C are similar, suggesting 
that the value of v has little effect on the long-term dynamics, although panel C shows slightly 
more variability for very small values of p > 0. 

5.2 Invariant densities about an attracting periodic solution 

Figs. [TTHm show two-dimensional invariant densities of Ni, N 2 and N 3 . By assuming ergodicity, 
invariant densities were computed on a 256 x 256 grid of x and y values from 10 8 consecutive 



Figure 9: Panel A shows the invariant density of N\ with the same parameter values as Fig. [8] 
and /i = 0.03. The value of the density is indicated by colour (dark red - the maximum value 
of the density; dark blue - zero). Panel B shows the invariant density at three times the noise 
amplitude as panel A. Panel C plots Oj - the fraction of instances that points return to x > 0 in 
j iterations (15. 2p - against the noise amplitude. 
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iterates of a single orbit with transient points omitted. 

Let us first explain panel C of Figs. I91HT1 Given a sample orbit of N\, N 2 or A 3 , for each point 
with Xi > 0 we let X, : be the smallest positive integer for which x l+ x, > 0, as in pTT] . T % represents 
the number of iterations required for a return to x > 0 from the point (xi,yi). Numerically we 
can compute a large number of values of (say M of them). Then for each j e Z + , we let <jj 
denote the fraction of the I, that are equal to j, i.e. 


aj ~ M Xj ~ Xi ’ 

i with. Xi >0 


(5.2) 


where Xk — 1 if A: = 0, and Xk — 0 otherwise. Figs. [911111 correspond to ft — 0.03 for which there 
is an underlying attracting 3-cycle. Therefore with small noise, <73 ~ 1, and for all j ^ 3, Oj ~ 0. 

Fig. [9] corresponds to the map N\. In panel A, the size of the noise is relatively small, so 
iterates of N± follow close to the 3-cycle. The invariant density is well-approximated by a scaled 
sum of three Gaussian densities centred at each point of the 3-cycle. About the point with 
x ~ —0.1, the density is stretched substantially more in x-direction than in the ^-direction. This 
is because points with x ~ — 0.1 have likely just undergone an iteration under the right half-map 
of Ni which is stochastic in the x-component but not the ^-component. The stretching around 
other iterates is then a consequence of iterating under N\ with x < 0 . 


A a = l B a = 3 



x x 



Figure 10: Panel A shows the invariant density of A 2 with the same parameter values as Fig. [ 8 ] 
and ji = 0.03. Panel B shows the invariant density with a = 3, and panel C is a plot of the 
fractions Uj (15.2[) . 
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With larger values of a, it is relatively common for the orbit to return to x > 0 in a number 
of iterations other than three, Fig. [9]-C. For this reason the invariant density displays additional 
characteristics. For instance with a = 3, the orbit returns to x > 0 in two iterations almost 10% 
of the time. Consequently, a substantial part of the invariant density centred roughly about the 
point of the 3-cycle with x ~ —0.04, lies in x > 0, Fig. [9]-B. The invariant density in panel B 
also has a small component with x & 0.05. This is due to points of the orbit with small values 
of x > 0 mapping under the left half-map of Ni due to the noise (i.e. returning to x > 0 in only 
one iteration). 

Fig. DU] illustrates A’ 2 using the same parameter values. Again with small noise the invariant 
density is roughly Gaussian about each point of the 3-cycle, whereas for relatively large noise 
iterates often cross into x > 0 prematurely causing the invariant density to take an irregular 
shape. When a = 3, points of the orbit that do not return to x > 0 in three iterations, usually 
return to x > 0 in two iterations. 

Fig. dU corresponds to N 3 and is similar to the previous figure. This indicates that the 
correlation time v has little effect on these pictures, although the invariant density has a slightly 
different shape when a = 3. 



Figure 11: Panel A shows the invariant density of N 3 with the same parameter values as Fig. [8] 
and ji = 0.03. Panel B shows the invariant density with a = 3, and panel C is a plot of the 
fractions (15.21) . 
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5.3 Invariant densities near the grazing bifurcation 

Fig. [HI shows invariant densities of Ah, N 2 and N 3 for parameter values closer to the grazing 
bifurcation than the previous three figures, specifically /i = 0.001. At this value of y, there is 
an underlying attracting 4-cycle, panel A. The noise amplitudes are given by (14. 7p . (14. lip and 
(I4.23P with a — 1. Recall, these amplitudes were chosen such that the size of the noise response 
of the three maps is roughly the same for larger values of y. Here, however, the size of noise 
response differs substantially. In panels C and D, which correspond to N 2 and N 3 respectively, 
the invariant density is approximately a scaled sum of four Gaussians about each point of the 
4-cycle. The invariant density in panel D, corresponding to v = 0, is noticeably larger than that 
of panel C. 

In panel B, which corresponds to Ah, the noise has a substantial effect because the switching 
condition of Ah is stochastic, and many points of Ah fall close to x — 0. The invariant density has 
a small C-shaped component in x > 0 corresponding to consecutive points of the orbit mapping 
under the left half-map of AR The part of the invariant density for x < 0 and y ~ 0.001 
corresponds to images of points under the right half-map of Ah, and is bimodal because the 
invariant density has roughly two components in x > 0. 



Figure 12: Panel A shows the attracting 4-cycle of (11.31) with the same parameter values as Fig. [2] 
and y = 0.001. Panels B, C and D show the invariant densities of Ah, N 2 and N 3 , respectively, 
using the same parameter values as Fig. [8] and y = 0.001. 
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5.4 Invariant densities about coexisting attractors 

With y = 0.0145 in Fig. [21 there is an attracting 4-cycle and an apparently chaotic attractor. 
These are shown in Fig. I~HHA. As with smooth maps [CT [52] . in the presence of noise orbits 
commonly dwell near the attractors for relatively long periods of time, and switch between 
attractors quickly. Invariant densities of N \, N 2 and N 3 are shown in panels B, C and D. In 
each case the bulk of the density is centred about the two underlying attractors. With white 
noise (panel D) there is no gap in the invariant density around (x,y) ~ (—0.12,0.012) due to 
randomness in both the x and ^-components of IV 3 . 

Lastly, Fig. [H] illustrates stochastic dynamics with y < 0. This figure corresponds to k osc = 5 
(different to Fig. [ 21 ) and y = —0.002 at which N has an attracting 3-cycle as well as the attracting 
fixed point ( x L ,y L ) (11. 5[) . These are shown in panel A. The dynamics of N 2 and N 3 for x < 0 
is deterministic, hence ( x L ,y L ) is a fixed point of these maps. Given an initial point (x 0 ,y 0 ) 
near the 3-cycle, sample orbits of N 2 and N 3 eventually reach ( x L ,y L ). In contrast, N { has an 
invariant density concentrated about the two attractors, panel B. The part of the density with 
x ~ 0 and y <yi corresponds to points of the orbit repeatedly following the left half-map of N\ 
(with Xi < 0 and Xi + ^3 < 0 ). 


A 

0.03 

y 

0.02 


0.01 





Figure 13: Panel A shows an attracting 4-cycle and a numerically computed attractor of (II. 3p 
with the same parameter values as Fig. [21 and y = 0.0145. Panels B,C and D show the invariant 
densities of N \, N 2 and IV 3 , respectively, with the same parameter values as Fig. [ 8 ]and y = 0.0145. 
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Figure 14: Panel A shows the attracting fixed point ( x L ,y L ) (II. 5p and the attracting 3-cycle of 
N with the same parameter values as Fig. [21 except k osc = 5 and fi = —0.002. Here r ~ 0.0927, 
5 ~ 0.1518 and x — 1- Panel B shows the invariant density of the corresponding map Ah, with 
v = 0.5 and e = 0.0005. 


6 Conclusions 

This paper concerns grazing bifurcations for which the associated dynamics is described by the 
Nordmark map (11.31) . The potential influence of randomness and uncertainties on the dynamics 
of (II.3j) was described in [37] by studying (ll.3j) in the presence of additive white Gaussian noise. 
Such a noise formulation is suitable if the nature of the randomness in the ODE system is 
practically independent to the state of the system, such as if there is a random forcing term. 

In this paper we considered the alternate scenario that randomness is present in the switching 
condition of the ODE system, and in /# - the part of the vector field opposite to the tangential 
intersection of the grazing periodic orbit. These cases are especially relevant for vibro-impacting 
systems for which impact events represent the primary source of uncertainty. We derived three 
different stochastic versions of (II.3p . These are the maps N\ (14.61) . which corresponds to a noisy 
switching condition in the ODE system, N 2 (14.101) . which corresponds to noise in /# with a large 
correlation time, and A r 3 (14.201) . which corresponds to white noise in f r. In each case the noise is 
nonlinear and non-additive. This indicates that some diligence should be taken when formulating 
stochastic return maps for grazing bifurcations of piecewise-smooth systems. 

The stochastic dynamics of Af, N 2 and N 3 differs in many ways to that of (11.31) with addi¬ 
tive noise, described in [ST]. For Ah, N 2 and AT 3 , dynamics prior to the grazing bifurcation is 
deterministic, and beyond the grazing bifurcation two-dimensional invariant densities are often 
skewed dramatically so that they appear almost one-dimensional. 

Near the grazing bifurcation, Af exhibits a large noise response relative to N 2 and N 3 . This 
suggests that if experimental data of a physical system shows relatively high variability near 
a grazing bifurcation, then it is likely to be most appropriate to include randomness in the 
switching condition of a mathematical model. Invariant densities of Af near grazing are highly 
irregular due to the randomness in the switching condition. For N 2 and N 3 , the size of the noise 
response increases, for most part, with the distance (in parameter space) beyond the grazing 
bifurcation. The maps N 2 and N 3 exhibit qualitatively similar invariant densities, which implies 
that the correlation time has little effect. Indeed the correlation time only influences the short- 
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time dynamics of (I4.8[) with u > 0, and the precise nature of these dynamics has a negligible 
effect on the invariant densities, which relate to long-time dynamics. 
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